Index theorem for topological heterostructure systems 
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We apply the Niemi-Semenoff index theorem to an s-wave superconductor junction system at- 
tached with a magnetic insulator on the surface of a three-dimensional topological insulator. We 
find that the total number of the Majorana zero energy bound states is governed not only by the 
gapless helical mode but also by the massive modes localized at the junction interface. The result 
implies that the topological protection for Majorana zero modes in class D heterostructure junctions 
may be broken down under a particular but realistic condition. 
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I. INTRODUCTION 



Zero energy bound states in vortex cores of super- 
conductors have been of much current interest in con- 
densed matter physics. Some classes of the vortex Ma- 
jorana states, obeying non-Abelian statistics, may serve 
as qubits for quantum computation.— — There are various 
theoretical proposals for realizing non-Abelian Majorana 
fermions in the core of vortices in topological supercon- 
ductors, e.g., a chiral p- wave superconductor, etc£~— Be- 
sides such vortex zero modes, topological superconductor 
junction systems, in which the order parameter changes 
sharply in real space, possess generically non-Abelian 
Majorana fermions^ A useful method for charactering 
the existence of the zero energy bound states localized 
at point defects such as vortices or point intersections 
consisting of the junction interfaces is the index theorem 
for an open infinite space, derived by Callias^i and by 
Weinberg^, and generalized by Nicmi-ScmenoffiiZ This 
theorem reveals the relationship between the zero energy 
bound states and the topology of background fields at 
large distance from the point defects. In this paper, we 
investigate the index theorem for the heterostructure sys- 
tem involving the topological insulator (TI), mainly fo- 
cusing on the supcrconductor-TI-fcrromagnct insulator 
junctions. We find that the number of Majorana zero 
modes is controlled not only by the phase winding of the 
superconducting gap, but also by non-topological mas- 
sive bound states localized at the junctions. 

The organization of this paper is as follows. In SeclJTl 
we first present our main results for the index of the 
superconductor-TI-ferromagnet insulator junctions, and 
discuss its physical implications. We give, in particu- 
lar, a physical explanation on how non-topological mas- 
sive bound states affect the index for Majorana zero 
energy modes. Our results are based on the cele- 
brated Niemi-Semenoff index theorem. To make this pa- 
per self-contained, we briefly review the Niemi-Semenoff 
index theorem in Sec lIIII In Sec lIVl we apply the 
Niemi-Semenoff index theorem to the superconductor- 
TI-ferromagnet insulator junctions, and obtain the in- 
dex theorem for topological heterostructure systems. In 



Sec[V] we also apply our results to topological insulator- 
ferromagnetic insulator heterostructure systems. We 
conclude in SeclVIIwith some discussions. 



II. SETUP AND MAIN RESULTS 

We consider the heterostructure system composed of 
an s-wave superconductor 7r-junction and ferromagnetic 
insulators placed on a TI, as depicted in Fig. [TJa), and in- 
vestigate the zero energy bound state localized at a point- 
like defect formed by the intersection of the 7r-junction 
interface and the ferromagnetic domain wall. The effec- 
tive Bogoliubov-de Gcnncs Hamiltonian is written as 



% = -ivTzo-jdj + AiTi + A 2 t 2 + h ■ a - /xt 3 , 



(1) 



where j = 1,2, r = (ti,t 2 ,t 3 ) and cr = (01, 0-2,03) are 
the Pauli matrices for the Nambu and spin space respec- 
tively, v is the velocity of the Dirac fcrmion, Ai and 
A2 are the real and imaginary parts of the gap function, 
h-cr is the Zeeman term, and /x is the chemical potential. 
Note that for the 7r-junction considered here, A2 = 0. 
It is also assumed that the thickness of the supercon- 
ducting film is sufficiently smaller than the penetration 
depth, and hence z-dependence of h is negligible. This 
system ([T]) belongs to class D in the Altland-Zirnbaucr 
symmetry classe s 18 ' 19 , and the vortex zero modes obeying 
non-Abelian statistics 9 - are classified as the Z2 invariant. 
The index theorem which is a main tool in this paper 
is applicable only to systems with chiral symmetry, i.e. 
UHlF = — TL is satisfied for a unitary operator IT. This 
symmetry is, however, not preserved for ([T}, because of 
/1303 and — /UT3 terms. Nevertheless, as will be clarified 
below, the Z2 invariant of ([T]) can be generically obtained 
from the index calculated for the case with ft, 3 = /j, = 0. 
Thus, we first neglect these two terms to calculate the 
index. We furthermore omit the h\U\ term to simplify 
the analysis, since this term does not affect the index 
of our system, as long as h\ is sufficiently small, and 
does not close the bulk energy gap. Then, the Hamilto- 
nian is reduces to that with chiral symmetry (class BDI), 



r 3 CT3'Hr30'3 = -ti, 



n 



-ivr^ajdj + AiTi + /12C2, 



(2) 



and its ground state is classified by Z. This enhanced 
topological number can be computed by the index the- 
orem described in detail in the next section. Note here 
that particle-hole symmetry, tio-{H,*T2<J2 = —T-L, valid 
for ([1]) as well as ([2]), ensures that the number of vortex 
zero modes are conserved modulo 2 even if the neglected 
chiral-symmetry-breaking terms are switched on again. 
Thus, the Z2 invariant of ([I} can be derived from the 
parity of the index of ([2]). 

Our central finding is that the heterostructure system 
composed of an s-wave superconductor 7r-junction and 
ferromagnetic insulators on a topological insulator (as 
shown in Fig (Ha)) described by ([2]) has the index: 



ind% =- [sign(/i+) - sign(/i_)] 



sign(h + )N x 



sign(/i_)7V x _ > _ c 



(3) 



where sgn(/i± ) is the sign of the asymptotic Zccman field 
hi{x — > ±00, y), and N x ^± oa are integer numbers which 
count how many times the band inversion occurs for mas- 
sive bound states at the 7r-junction, as the Zeeman mag- 
netic field increases from zero to h-2(x — > ±00, y). (see 
Fig. [2] and discussion given at the end of this section.) 
Note that in Eq.Q, we take the origin of the xy co- 
ordinate (x, y) = (0,0) at the location of the point de- 
fect in FigJTJa). It is also naturally assumed that the 
sign of h,2(x — > ±00, y) is independent of y. Especially 
in the cases of the uniform asymptotic Zeeman field, 
h(x —> ±00, y) = h±, the index (|3|) is simplified to 



ind-H = - [sign(/i+) - sign(/i_)] 



sign(>+) Yl - si S n (^-) 5Z 

E n <\h+\ E n <\h-\ 



(4) 




m<0 
Topological Insulator 



I/!- 



FIG. 1: (a) The heterostructure geometry for an s-wave su- 
perconductor (SC) 7r-junction and a ferromagnetic insulator 
(FMI) on the surface of a topological insulator. A filled cir- 
cle at the interface is a point defect formed by the intersec- 
tion of the 7r-junction and the ferromagnetic domain wall, 
(b) The heterostructure geometry for a topological insulator- 
ferromagnetic insulator tri-junction. 



where h+ (/i-) is a Zeeman field at the 7r-junction in- 
terface for x > (x < 0), which induces mass gap of 
the one-dimensional gapless helical Majorana mode lo- 
calized at the junction interface, and E n (> 0) denote the 
absolute value of the mass gaps of the one-dimensional 
massive modes localized at the junction interface. The 
sum in ^ is taken only for one part of the Kramer's pair. 
As mentioned before, the sum in ^ represents the num- 
ber of times the band inversion occurs for the massive 
bound states, as h± increase from zero to finite values. 
Because of particle-hole symmetry, this counting can be 
expressed only by E n > 0, as shown in Eq. (|1]). 

The index ((4]) (or ([3])) expresses the number of zero 
energy Majorana bound states in a point-like defect at 
the junction of a chiral-symmetric superconductor (class 
BDI in the Altland-Zirnbaucr symmetry classe a 18 ' 19 ). 

The index ([4J is interpreted as the phase winding of 
the superconducting gap A around the point defect. In 
the case of the 7r-junction with a Zeeman field as shown 
in the FIG. QJa), the change of the phase of the gap 
function can be defined in the following way by using the 
Teo-Kane's adiabatic argument i 21 ' 22 
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FIG. 2: (Color online) A schematic picture of the band in- 
version of massive bound states, (a) Black lines denote the 
energy band of helical Majorana fermion, while red and blue 
lines denotes that of massive bound states. Lower panels indi- 
cate the right-half part of the heterostructure geometry shown 
in FIG. [Ha). 

Without a Zeeman field, the 7r-junction possesses a he- 
lical Majorana fermion localized at the junction interface 
(FIG. Eta))j*i The Zeeman field from the ferromagnetic 
insulator lifts the Kramer's degeneracy, and induces a 
mass gap of the helical Majorana fermion (FIG. HJb)). 
Adiabatic deformation of the Hamiltonian without clos- 
ing the energy gap enables us to introduce a nonzero 
imaginary part of the superconducting gap A2 at the 
junction interface. In this process, the sign of A2 is deter- 
mined by a Zeeman field such that sign(A2) = — sign(/i2). 
Hence, the phase shift is — n sign(/i2), which is described 
by the first term in Q. This contribution depends only 



on the sign of Zccman field and docs not depend on the 
detail of the junction interface. 

The second term in (0| is, on the other hand, a new 
contribution which was not discussed in previous litera- 
tures in the context of hctcrostructure systems and de- 
pends on the detail of the junction interface. The super- 
conducting gap Ai changes its sign at the 7r-j unction. If 
the spatial variation of the magnitude of Ai in the vicin- 
ity of the junction is sufficiently slow, there exist massive 
bound states localized at the junction, which come in 
Kramer's pairs with a mass gap \E\\ (FIG. [UJa)). The 
Zeeman field parallel to y-axis shifts the mass gaps of 
the Kramer's pairs by \E\ ±/i2|, respectively(FIG.[2jb)). 
(Here, we assume the Zeeman field is uniform. If not, the 
mass gap of the bound states depends on the detail of 
Zeeman field. But the qualitative nature is unaffected.) 
When I/12I reaches |-Ei|, the energy gap at the junction 
interface is closed (FIG. [U[c)), and a band inversion oc- 
curs for /i2 < — \E\\, After band inversion, the junc- 
tion interface structure acquires the — 2-k sign(/i 2 ) phase 
shift in addition to the — it sign (^2) phase shift, result- 
ing in the total — 3w sign(ft_2) phase shift (FIG. HJd)). 
This additional 2-k phase production arises for each mas- 
sive bound state with mass gap \E 2 \, I-E3I, ■ ■ ■ , which de- 
scribes the second term in Q . Therefore, massive bound 
states at the interface give rise to additional phase wind- 
ing around the point defect formed by the intersection 
of the 7r-junction interface and the ferromagnetic domain 
wall shown in FIG. [lja). 

This new contribution from the second term of (U]) 
has an important implication for the class D hetcrostruc- 
ture system. As mentioned before, the class D system is 
characterized by the Z2 invariant for the Majorana zero 
modes, which is exactly the parity of the index Q ob- 
tained by switching off chiral-symmetry-breaking terms. 
Thus, the Z2 invariant of the class D hctcrostructure 
system may be changed by the non-topological massive 
bound states. This leads to a breakdown of topological 
protection of Majorana vortex modes when the second 
term of (@| is an odd integer. 



III. NIEMI-SEMENOFF INDEX THEOREM 



with V = —ijidi+K(x), for the basis that the T-matrices 
are represented as 



" ' °I 0.1 



1 
-1 



(0) 



Here, the indices i = l,2,...,d are those for the spatial 
coordinates, the % matrices are constant matrices that 
satisfy 7,7] + 7,- 7? = 7] 7,- + 7)7, = 2S tJ , and K (x) in- 
cludes all background fields such as electromagnetic fields 
and the superconducting gap. The index of the Hamil- 
tonian is defined by ind'H := dim ker 2?' — dim ker X>, 
which is the difference between the number of zero en- 
ergy states of H. with the opposite chirality. We assume 
all background fields are asymptotically independent of 
the normal coordinate, hi(x)diQ(x) — s> (|a;| — > 00), 
where h{x) is a unit vector normal to an open boundary 
at J x J — »■ 00. It is known that ind T-L is expressed as the 
sum of the volume integral of the chiral anomaly and the 
surface integral of the chiral current , 16 ' 17 



ind"W 

= / d d x tr (x I T 5 I x) + 



dS tr 



it{x)T 5 H- 



(7) 



where dS is the volume element of the boundary, and 
T(x) := fii{x)Ti. The definition of the terms of r.h.s. 
in ([7]) needs appropriate regularization. In this paper, 
we symbolically use the expression of r.h.s. in ([7]). The 
first term in ([7]) is the integrated chiral anomaly which 
is present only in even spatial dimensions. When d = 2, 
it is explicitly written in terms of the background field 
Q(x) asi 6 - 



d z x tr (x J T 5 I x) = -^- I d z xtr iY 5 T i d l Q{x). (8) 



This formula will be used later, (see Eq.tfTB"]) below) The 
second term in (jTj) is the boundary integral of the chiral 
current density normal to the boundary, and this term 
can be rewritten as the spectral asymmetry constructed 
from the real part of the eigenvalues of a certain bound- 
ary operator M. as shown below, 



In this section, for the convenience of readers, we 
briefly review the Niemi-Semenoff index theorem which 
is used for the derivation of our results in the following 
sections. The Niemi-Semenoff index theorem relates the 
number of zero energy modes in Dirac fermion systems to 
the geometrical structure of spatially varying mass terms. 
In particular, the index is determined by the asymptotic 
behaviors of mass terms at open boundaries. We con- 
sider the Dirac Hamiltonian with chiral symmetry in d- 
dimensional space with open boundaries, the Hamilto- 
nian of which is given by, 



U = -iTidi + Q(x) = 







(5) 



tr (x ir(x)r 5 ^" 






tr ( x 



tr ( x 



*/(x) 
yt(x) 



(i^(x)V) 



1 
-1 



(i^(x)V) 



V 

pt 



(9) 



Here 7(x) = fn{x)^i are the normal components of 7 ma- 
trices, and we used the cyclicity of trace. The boundary 
operator M. is defined by irf {x)T> = d + Ai , 



M = ^{x)^Jd % + irf(x)K(x), 



(10) 



where jfix) = 7i — j(x)fii(x) are the tangential compo- 
nents of 7 matrices, and d = hi(x)di is the directional 
derivative normal to the boundary. We assume % does 
not possess zero modes at infinity, which corresponds to 
the absence of zero modes in Ai. Since Ai is indepen- 
dent of the coordinate normal to the boundary, we can 
introduce the Fourier transformation for the normal co- 
ordinate: 



- (p dStr ( x 



iTT 5 H- 



- (j> dSti < x 

— [ dki 
4tt ./_«, J 



1 



1 



M + d M^-d 
dStr ( x 



1 



M+ik M^—ik 



(11) 



Introducing the cigenmodes Aitfi — X(j) and M'i/> = A*?/;, 
we rewrite Eq. (|TT|) as 



1 
-in 



dX p(A)sign [Re(A)] 



1 



A* — ik 



(12) 



- v (Bb(M)). 



where p(A) is the spectral density of boundary operator 
Ai. This term is the spectral asymmetry constructed 
from the real part of the eigenvalues of M.. Eventually, 
ind % is written as^£ 



This is the Niemi-SemenofF index theorem for an open in- 
finite spaced The integrand of the anomaly contribution 
is generally the total derivative. Hence ind H. depends 
solely on the asymptotic behavior of background fields. 



IV. MAJORANA ZERO MODES AT A POINT 

DEFECT IN 

SUPERCONDUCTOR-FERROMAGNET 

INSULATOR HETEROSTRUCTURE SYSTEMS 



In this section, we derive the index ^ for the topo- 
logical heterostructure system depicted in Fig. QJa) by 
applying the Niemi-Scmcnoff index theorem explained in 
the previous section. For this purpose, we first obtain the 
boundary operator (fTTJj) for the Hamiltonian ©. This is 
achieved by the following procedure. By applying the 
unitary transformation, 



with 



T3C3 



a 3 
-0-3 



H^ 



Ur 3 a s U^ 



1 
-1 



(14) 



U 



(15) 



indU = / d d x tr (x\T 5 \x) + -r\ (Re (M)) . (13) the Hamiltonian © is represented as 



H^ 



-v(a 2 d x 





Oldy) 



h 2 (J2 + Ai 



v(a 2 d x 



endy) 




h 2 <T 2 + Ai 



(16) 



r 



In this representation, (71,72) = (ic2, —i&i), and 
K(x,y) = h 2 a 2 + Ai. Then, the boundary operator .A/f 
defined by JTUJ) is 



M = wa^dx + Aictt + hr ~ iha 3 



(17) 



where a = niai and ax = nfa.i are components of a 
vector a = (ai,a 2 ) which are, respectively, normal and 
tangential to an open boundary at |a;| — > 00 or |y| — ^ 00. 
(We note that the origin of the xy coordinate (x, y) = 
(0, 0) is taken at the position of the point defect in Fig. 
[IJa).) From ©, the anomaly part of the index is 



/ d 2 xtr (x I T 5 I x) = - Id' 



X^ijOifi j 

7T 



dlh. 



(18) 



To simplify the analysis, we consider the kink structure 
for the gap function Ai(y) = Atanh(y/£) with A > at 
the interface of 7r-junction. As we shall see momentarily, 
the parameter £ in the gap function describing the width 
of the kink gives rise to a crucial effect on the index. 
Here, it is natural to assume that the asymptotic value 
of the magnetic field, h 2 (x,y) — > h±(y) as x — > ±00, 
have nonzero values with definite signs, sign(h±), near 
the interface of the 7r-junction. On the other hand, the 
magnitudes of h±(y) at y — > ±00 does not affect the 
index of our system. Thus, for simplicity, we assume 
h±{y) — >• for y — >• ±00. 

Now we calculate the spectral asymmetry 77(Rc(A^)) 
for the heterostructure system depicted in Fig. QJa). It 
follows from (|17|) that the boundary operators Ai for 



±00 and y — > ±00 arc given by, 



M x . 
M x . 

My- 



Xv) = ~iva 3 d y + Ai(y)cr 2 + h + (y), 
oo(y) = -ivcr 3 dy - A 1 (y)a 2 - h-(y), 






-iva 3 d x 
iva 3 d x 



-Aa u 
Aa ± . 



(19) 
(20) 

(21) 
(22) 



The spectral asymmetry (|12[) is the sum of the par- 
tial spectral asymmetries of four sides. Each bound- 
ary operator is hermitian, since the Zeeman field normal 
to the boundary in (JTTJ) vanishes. The spectral asym- 
metries for A4 y ^±oo are zero, since the eigenvalues of 
■A4j,->±oo come in pairs ±A due to "chiral" symmetry 
02 A^-*. ±00 02 = — 7W a _j.±oo- To calculate the spectral 
asymmetries for Adx-*±oc, we exploit an approach devel- 
oped by Lott : 23 ' 24 Let H T be a one-parameter family of 
Hamiltonians defined in r G [0, 1] which interpolate be- 
tween a reference Hamiltonian Ho and Hi = M x ^±oo- 
In the calculation of the spectral asymmetry, we choose 
the reference Hamiltonian Ho for which the spectral 
asymmetry r\ ("Ho) is known. The variance of the spectral 
asymmetry as a function of r is composed of two parts: 
one is the continuum part vfc raised by the change of 
high energy continuum energy spectrum, and the other 
is a discrete part arising from the spectral flow which 
changes by ±2 when a discrete eigenvalue A„(r) crosses 
zero from negative (positive) to positive (negative) en- 
ergies: A [sign(A„(r))] = ±2. Thus, we can write the 
spectral asymmetry in the form, 



r,(H) = 7,(7*0) 



dn c 
dr-^ + 2 (spectral flow). (23) 

(XT 



Let us consider the spectral asymmetry for x — > 00. 
In this case, the boundary operator A^-kx) is basically 
the Jackiw-Rebbi Hamiltonian. Therefore, in the ab- 
sence of a magnetic field, h + (y) = 0, M x ->oo possesses 
the Jackiw-Rebbi zero mode localized at the interface 
of 7r-junction, <j> Q {y) ex '(1, -l)e~ /" **(v')Mv' . When 
a magnetic field is switched on, a finite value of h+(y) 
shifts the bound state energy from zero to a nonzero 
value with the same sign as h + . The spectral asym- 
metry for this boundary operator was previously com- 
puted by Lott^ 3 -. Using the reference Hamiltonian given 
by Ho = iva 3 d y + Ai(y)a 2 — osign(/i+)0i with a small 
positive constant 8, we obtain, 



1 



v(M x ^oo) 



1 If 

-sign(/i + ) / dyh + (y) + (spectral flow). 

2 it J 



(24) 



The first term conforms to the fermion fractionalization 
in the Jackiw-Rebbi system . 25 ' 26 The second term is the 
volume part of the variation of r\ T , which cancels out the 
anomaly contribution of the index (fT8|) . The third term is 
the spectral flow contribution from H to .M :r _s. 00 which 
depends on the structure of Ai(y) and h + (y). The spec- 
tral flow stems from the bound states localized at the 




Continuum states 



FIG. 3: The r-dependent energy spectrum of a one parameter 
family of Hamiltonians H. T : H T — (1 — t)Ho + tH. This figure 
shows the case of £ = 5£ c , h+ = 0.7A > 0, and 5 = 0.05A. 



interface of the 7r-junction. Since the Lott's derivation 
of (|24|) in ref[23| is highly technical, we give a more ele- 
mentary derivation of (|24[) in Appendix, which we believe 
is useful for readers. If the finite value region of ft+(y) 
is much longer than £, h+(y) is approximated as a con- 
stant chemical potential, h + (y) — > h + . In this case, the 
y-dependent part of A^-xjo, i.e., iva 3 d y + Atanh(y/£) is 
exactly solvabl e 27 ' 28 The eigenvalues of this Hamiltonian 
are E = 0, 



E 



n,± 



±A 



1,2,.. 



,<!/), (25) 



and E p< ± = ±y / v 2 p 2 + A 2 , (p e R), where v = £A/v = 
£/£ c is a ratio of £ to the coherence length of the su- 
perconducting state : £ c = v/A. Eq and E n> ± are the 
energy of the bound states localized at the interface of 
7r-junction. The Jackiw-Rebbi zero energy bound state 
exists for arbitrary £, while the massive bound states ex- 
ist only when £ > £ c . h + induces constant shifts to eigen- 
values (|25|) . and hence, the spectral flow is given by the 
bound states that cross zero between E n ± and E n ± + h + 
as shown in Fig. [3] 



(spectral flow) 



sign(/i+) Yl ' 

E n <\h+\ 



(26) 



where E n = E n ,+ - Note that due to the term 
— (5sign(/i + )(Ti in the reference Hamiltonian Ho, the spec- 
tral flow from the Jackiw-Rebbi zero energy bound state 
is excluded in the sum (|26[) . In a similar way, the spectral 
asymmetry of At^^-oo, ^(Adx-^-oo) can be calculated 
as 



j 1 /" 

--sign(/i_) + - / dyh_(y) - sign(/i_) ^ 



(27) 



E„<\h- 



Using Eqs. CCD, dHJ), ([Ml), flU and (J^TJ) together, we 
arrive at the formula ((4]). 

So far we have calculated the index for $Z§. The 
Z2 index N for the vortex zero modes in class D het- 
crostructurc system (TT]) is, as we have mentioned, given 



by N = indH (mod 2). Remarkably, the second term 
in @, which is basically the contribution from non- 
topological bound states of the Jackiw-Rcbbi Hamilto- 
nian, can affect the Z 2 index. This indicates that the ex- 
istence of the non-Abelian vortex zero modes depends on 
the kink-structure of the gap function parametrized by £. 
Actually, in the case that the second term of (fj| is equal 
to an odd integer, which can indeed occur when £ > £ c 
(i.e. v > 1), the Z 2 invariant for Majorana zero modes 
is changed, resulting in the breakdown of the topological 
protection for the Z 2 Majorana modes. We now discuss 
the condition for which £ > £ c is realized. Actually, to de- 
termine £ precisely, we need to solve the Bogoliubov-de- 
Gennes equation for proper boundary conditions, which 
is out of the scope of this paper. Instead of present- 
ing such precise analysis, we here give a qualitative ar- 
gument. For the superconductor-ferromagnetic insulator 
junction as depicted in Fig. 1 (a), the gap function at the 
junction is reduced by magnetic scattering at the inter- 
face between the superconductor and the ferromagnetics 
On the other hand, the dimension of the ferromagnet 
insulator along the y-axis denoted as L y plays the role 
of the characteristic length scale for the spatial varia- 
tion of the exchange field along the y-axis. Thus, when 
L y is sufficiently larger than the coherence length £ c , we 
can neglect the spatial variation of the superconducting 
gap raised by magnetic scattering near the interface, and 
hence it is expected that £ < £ c is satisfied, ensuring 
the topological protection of Majorana modes. However, 
when L y is comparable to £ c , the spatially inhomoge- 
neous reduction of the superconducting gap due to mag- 
netic scattering crucially affects the magnitude of the pa- 
rameter £. In particular, when L y is slightly larger than 
£, it may be possible that £ > £ c is realized, which leads 
to the above-mentioned mechanism of the breakdown of 
Z 2 nontriviality. 



V. ZERO MODES IN LINE DEFECTS OF 

TOPOLOGICAL INSULATOR-FERROMAGNET 

INSULATOR HETEROSTRUCTURE JUNCTIONS 

The index theorem Q is also applicable to a topologi- 
cal insulator-ferromagnetic insulator tri-junction system, 
the setup of which is depicted in Fig. [1] (b). The two 
orbital effective Dirac model for this system is written 
as2C 

H = vk z [i 1 o- i - ivuxdjdj + m{y)iJ,3 + h 2 {x 1 y)a 2 - n, 

(28) 

where j = 1,2, (J, = (/ii,/x 2 ,/x 3 ) and <x = {a u a 2 , 0-3) are 
the Pauli matrices for the orbital and the spin spaces, 
respectively, v is the velocity of the Dirac fermion, m is 
the mass gap whose sign determines whether the sys- 
tem is in a topological (m < 0) or trivial (m > 0) 
phase, h ■ er is a Zccman term, and [i is the chemical 
potential. We have assumed the translational invariance 
along the z-direction. This system belongs to the class A, 



and the chiral gaplcss modes localized at line defects are 
classified as Z3i As in the case of the superconductor- 
ferromagnet insulator junction, we, first, neglect chiral- 
symmetry breaking terms, putting /i = 0: 



H = vk z ^a z - ivni<Tidi + m(y)/x 3 + h 2 (x,y)a 2 
=: vkzHia 3 +H(x,y), 



(29) 



Because of chiral symmetry {/11CT3, T-L{x, y)} = 0, the chi- 
ral zero bound states of it{x,y) with chirality ± corre- 
spond to the chiral gapless modes with the energy dis- 
persion ±wfc 2 with chirality ±. Since H(x,y) is of the 
same form as @, indH is given by 0, but in this case, 
E n (> 0) is the mass gap of the two-dimensional massive 
bound states localized at the surface of topological insu- 
lator. The first contribution in (QJ agrees with the wind- 
ing of the Axion vortexiCC The second contribution in ^ 
correspond to the non-topological integer part of Axion 
field which depend on the microscopic structure of the 
interface between the topological insulator and the triv- 
ial insulator. The chiral gaplcss mode cannot be massive 
since backward scatterings are suppressed. This index 
for chiral gaplcss modes in class A survives against any 
perturbations. 



VI. CONCLUSION AND DISCUSSIONS 

Here, we remark a topological property for the in- 
dex of heterostructure systems. The index Q is stable 
against the continuous change of Hamiltonian % unless 
the boundary operator M. have zero modes, i.e., unless 
the Hamiltonian H have no gapless modes at infinity. 
In this sense, the index (g]) is topologically protected by 
the energy gap at the boundary. This feature is similar 
to the topological order in bulk systems protected by a 
bulk energy gap. 

In conclusion, we have shown that the number of Ma- 
jorana bound states in the 7r-junction-fcrromagnet het- 
erostructure systems is affected by massive bound state 
localized at the interface, which has an important impli- 
cation for topological protection of zero modes in class D 
systems. 
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Appendix A: Derivation of (|24[) 



In this appendix, we derive Eq. (|24p by using 
the Nicmi-Scmcnoff formula for the fermion number 
fractionalizationf24 which is an alternative expression of 
the Nicmi-Scmcnoff index theorem. 



1. Niemi-Semenoff formula of fermion number 
fractionalization 



2. The spectral asymmetry of the Hamiltonian (|19[) 
and (T20J) 



by 



We introduce an extended Hamiltonian "H(r, x) defined 



H(T,x) = -iY j T i d i + Q(T,x) 



o d T + ff (t, x; 

-d T + H{T,x) 



(Al) 



where <9o = <9 T , and H(r,x) = —iJ2i=ili^i + ^'( t j x ) 
satisfies H(t — > oo, x) = Hi(x), and -ff(r — > — oo,x) = 
Hq(x) with Hi(x) the target Hamiltonian for which we 
want to calculate the spectral asymmetry, and Hq(x) 
a reference Hamiltonian. -ff (r, a;) interpolates between 
Hi(x) and Hq(x) as a function of the auxiliary parame- 
ter r. Also, a new gamma matrices 70 = i is introduced. 
From the Niemi-Semenoff index theorem (fT3|) . the index 
of % is given by 

ind % = / drd d x tr (t, x \ T 5 \ t, x) + -77 (Re(X)) . 

(A2) 

Note that there are three boundary operators which 
contribute to the spectral asymmetry ?/(Re(A / f)); i.e. 
A4 T ->oo; A / ( r -s.-oo, and A"! i^^oo. The first two boundary 
operators arc, respectively, related to Hi(x) and Hq(x), 

Mr^oo = tfudi + itfK(r ->■ oc,x) = i?i(x), 

Mt^-oo = tffidi + i^K(r ->■ -00, x) = -H (x), 

(A3) 

where, 7(r — > ±00) = ±70- Then, from Eq. ([A2[) . we 
have^i 

- / rfrd d x tr(r,x|r 5 |r,x) - -77 (Re(X| a! |_ > . 00 )) , 

(A4) 

where we have used r)(—H) = —r](H). Generally, the 
change of the spectral asymmetry is divided into its con- 
tinuous part and discontinuous part, 

r}{H x ) - r](H ) = I dr^L + 2(spectral flow). (A5) 

Comparing (|A4[) and (|A5|) . we find, 

(spectral flow) = ind H, (A6) 

Then the spectral asymmetry of H\ is given by 

-77 (Hi) = -77 (H ) + (spectral flow) 

- / drd d x tr(r,x|r 5 |r,x) - -77 (Re(M\ a> \- > . 00 )) . 

(A7) 



We apply (|A7|) to the boundary operator (|T9|) . 

Hi(y) = iva 3 d y + A t (y)a 2 + h+(y), (A8) 

and the reference Hamiltonian, 

H (y) = iv<r 3 d y + Ai(y)cr 2 - Sai, (A9) 

where 6 is a small constant which is introduced to sup- 
press zero energy modes of Ho(y). We introduce a Hamil- 
tonian that interpolates Hi(y) and Ho(y), 

H(r,y) = iva 3 d y + Ai(y)a 2 + A 2 (r, y)cri + h 2 (r,y), 

(A10) 

where A 2 (t — > 00, y) = 0, A 2 (t — > —00, y) = —5, h 2 (r — > 
00, y) = h+(y), and h 2 (r —> —oo,y) = 0. We assume 
A(r, y) and h 2 (r,y) form a single kink structure along 
the r direction. The boundary operator at y — > 00 is 

■My— >-oo 

= 7*74^ + i^K(t, y -> 00) 
= -ia 3 d T - ia 3 K(r, y ->■ 00) 

= -ia 3 d T - Ai(cxd)(Ti + A 2 (r, oo)cr 2 - i/^r, 00)0-3, 

(All) 

where 7(3/ — > 00) = 71 = —(73. The fourth term in 
(jAlip does not contribute to the spectral asymmetry, 
since —ih 2 (r,oo)a 3 is anti-hcrmitc. Thus, the spectral 
asymmetry i] (Re (M. y -^oo)) arises from the Jackiw-Rcbbi 
Hamiltonian, — ia 3 d T — Ai(oo)<7i+A2(t, oo)cr 2 . The spec- 
tral asymmetry of this Hamiltonian is well known,— and 
equal to the phase winding of — Ai(oo)+7A 2 (r, 00) raised 
by changing r from —00 to 00 : 



-77 (Re (My-**,)) = -—Arctan 

Z Z7T 



(A12) 



, < Arctan < -|. 



,Ai(oo) 

where Arctan has principal values, - 

Similarly, the spectral asymmetry of the boundary oper 

ator at y — > — 00 is computed as 

1 



-rj(Re(M y ^- 00 )) 



_L Alcta „ ( « 

2?r v a i(-°°; 



(A13) 



On the other hand, the spectral asymmetry of the refer- 
ence Hamiltonian H is 



h™-h 



Arctan 



Ai(oo) 



Arctan 



Ai(-oo) 



S ) \ 5 

(A14) 

Summing up the boundary contributions (|A12|) . (|A13j) 
and the contribution from the reference Hamiltonian 
(| A14|) . we obtain, 



1 



1 



77 (H ) - -77 (Re(M y ^oo)) - -77 (Re(M y ^oo)) 



sgn 



Ai(oo) 



sgn 



2 

Ai(-oo) 



(A15) 



where we have used Arctan(2:)+Arctan(.T~ 1 ) = : |sgn(a;). 
If Ai(oo) > 0, Ai(-oo) < 0, and sgn(<5) = sgn(/i+) 
(sgn(/i+) is the sign of the Zeeman field h+(y)), then 
(|A15|) is written as, 



Prom (|A"7|) . (jA"l"6l) and (|ATf|) . we obtain (f24|) . 



-r? (flb) - ^ (Re(M^oo)) - ^ (Re(A^ i 



-)) 



-sgn(ft,_| 



(A16) 



Next, the anomaly part in (|A7|) is calculated from ©, 
/ drdy tr (r, y \ T 5 | r, y) = -— J drdytr [iT b TidiQ{T, y)) 
= -^JdTdy(-4d T h 2 (T,y)) 
= ~ I dyh + (y). 
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